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ABSTRACT 

We describe the discovery of the longest microlensing event ever observed, 
OGLE-1999-BUL-32, also independently identified by the MACHO collab- 
oration as MACHO-99-BLG-22. This unique event has an Einstein radius 
crossing time of 641 days. The high quality data obtained with difference im- 
age analysis shows a small but significant parallax signature. This parallax 
effect allows one to determine the Einstein radius projected onto the observer 
plane as te « 29.2AU. The transverse velocity projected onto the observer 
plane is about 79kms _1 . We argue that the lens is likely to be have a mass of 
at least a few solar masses, i.e., it could be a stellar black hole. The black hole 
hypothesis can be tested using the astromctric microlensing signature with 
the soon-to-be installed Advanced Camera for Surveys on board the Hubble 
Space Telescope. Deep X-ray and radio images may also be useful for revealing 
the nature of the object. 
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1 INTRODUCTION 

Gravitational microlensing is rapidly becoming an important astrophysical tool (for a review, 
see Paczynski 1996). The unique strength of this technique is that it provides a mass-selected 
sample for a variety of astrophysical applications, such as studying the Galactic structure 
and mass functions in the local group. So far, over one thousand microlensing events, mostly 
toward the Galactic bulge, have been discovered (e.g. Alcock et al. 2000; Wozniak et al. 
2001; Bond et al. 2001). Most (~ 90 per cent) microlensing events are well described by the 
standard shape (e.g. Paczynski 1986). Unfortunately, from these light curves, one can only 
derive a single physical constraint, namely the Einstein radius crossing time, which involves 
the lens mass, various distance measures and relative velocity (see §3). This degeneracy 
means that the lens properties cannot be uniquely inferred, thus making the interpretation 
of the microlensing results ambiguous. 

The parallax microlensing events are one class of exotic microlensing events that allow 
this degeneracy to be partially removed. The parallax effect we discuss here arises when the 
event lasts long enough that the Earth's motion can no longer be approximated as rectilinear 
during the event (Gould 1992; see also Refsdal 1966 for a related effect). Unlike the light 
curves for the standard events which are symmetric, these parallax events often exhibit 
asymmetries in their light curves due to the motion of the Earth around the Sun. These 
events allow one to derive the physical dimension of the Einstein radius projected onto the 
observer plane and hence the lens degeneracy can be partially lifted. 

A number of parallax microlensing events have been reported in the literature (Alcock 
et al. 1995; Mao 1999; Soszyhski et al. 2001; Bond et al. 2001; see also Bennett et al. 1997). 
Smith, Mao & Wozniak (2001) recently developed a method to systematically search for 
parallax signatures in the OGLE-II microlensing candidates found by Wozniak et al. (2001). 
We have uncovered several parallax candidates in this database. One of these, OGLE-1999- 
BUL-32, turns out to be the longest microlensing event ever observed. The purpose of this 
paper is to analyze this unique event in some detail. We argue that this event is likely to 
be caused by a stellar mass black hole; other black hole candidates from microlensing have 
been reported in conference abstracts (Bennett et al. 1999; Quinn et al. 1999). The outline 
of the paper is as follows. In §2 we briefly describe the observations, data reduction and 
our parallax search algorithm, in §3 we describe our model for this spectacular microlensing 
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event, and in §4 we propose future observations that can further test our model, particularly 
with the Hubble Space Telescope (HST). 

2 OBSERVATIONS, DATA REDUCTION AND SELECTION PROCEDURE 

All observations presented in this paper were carried out during the second phase of the 
OGLE experiment with the 1.3 m Warsaw telescope at the Las Campanas Observatory, 
Chile. The observatory is operated by the Carnegie Institution of Washington. The telescope 
was equipped with the 'first generation' camera with a SITe 3 2048x2048 pixel CCD detector 
working in the drift-scan mode. The pixel size was 24/xm, giving the scale of 0.417"per pixel. 
Observations of the Galactic bulge fields were performed in the 'medium' speed reading 
mode with the gain 7.1 e - ADU -1 and readout noise about 6.3 e~. Details of the instru- 
mentation setup can be found in Udalski, Kubiak & Szymahski (1997). The majority of the 
OGLE-II frames were taken in the /-band, roughly 200-300 frames per field during observing 
seasons 1997-1999. Udalski et al. (2000) gives full details of the standard OGLE observing 
techniques, and the DoPhot photometry (Schechter, Mateo & Saha 1993) is available from 
the OGLE web site at http://www.astrouw.edu.pI/~ ogle/ogle2/ews/ews. html. 

Wozniak et al. (2001) searched for microlensing events in the three year OGLE-II bulge 
data analyzed using difference image analysis. The difference image analysis pipeline is de- 
signed and tuned for the OGLE bulge data (Wozniak 2001), and is based on the algorithm 
from Alard & Lupton (1998) and Alard (2000). The difference image analysis pipeline re- 
turned a catalog of over 200,000 candidate variable objects, from which 520 microlensing 
candidates were identified using a combination of an algorithmic search, visual inspections 
and a cross-correlation with the candidates identified by Udalski et al. (2000) from the 
DoPhot analysis. The details can be found in Wozniak et al. (2001) and will not be repeated 
here. 

We then searched for parallax microlensing events using the method developed in Smith 
et al. (2001). Here we outline the prescription. In the first step, we fit each microlensing 
light curve with both the standard model and the parallax model (see §3 for the procedure 
applied to OGLE-1999-BUL-32). The events that show significant improvements with the 
incorporation of the parallax effect are then recorded and subjected to further studies. 
Among the recorded events, we then select those events for which the peak is at least 30 times 
higher than the noise level and the time interval during which the microlensing variability is 
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at least 3c above the noise level is longer than 100 days. These two filters properly account 
for the fact that (subtle) parallax signatures are most likely to be detectable in long-duration 
events and those with high signal-to-noise ratio. We found this prescription to be successful. 
Several good candidates and a number of marginal ones were identified. We refer the readers 
to Smith et al. (2001) for further details. 

In this algorithmic search, one microlensing event in Wozniak et al.'s catalog, sc33_3764, 
passed all our criteria. The microlensing variability was in fact first identified in real-time 
by the MACHO alert system as MACHO-99-BLG-22; it was also detected by the OGLE 
early-warning-system as OGLE-1999-BUL-32. The star, however, first escaped detection as 
a microlensing candidate (or even as a 'transient', see Wozniak et al. 2001) in the differ- 
ence image analysis, because the star never reached a 'constant' baseline (see Fig. 2). The 
event was recovered by cross-correlating the variable stars with the microlensing candi- 
dates found by Udalski et al. (2000). Throughout this paper, we shall refer to this event 
as OGLE-1999-BUL-32, following the notation of Udalski et al. (2000). The position of 
the star is RA=18:05:05.35, and DEC=-28:34:42.5 (J2000). The Galactic coordinates are 
/ = 2°. 460, b = —3°. 505. The DoPhot photometry and finding chart for the star are available 
online, t The total /-band magnitude of the lensed star and the nearby blend(s) is about 
/ pa 18.1 (uncertain by about 0.05mag, Wozniak et al. 2001). The baseline magnitude of the 
lensed star alone is about I ~ 19.2 (see §3). There are several U-band frames in the OGLE- 
II database when the composite was fainter than / = 16.6 magnitude. The average V — I 
colour of the composite is about 1.6. Fig. 1 shows the colour-magnitude diagram for the stars 
within a field of view 3.8' x 8' around OGLE-1999-BUL-32. From this figure, it is clear that 
the magnitude and colour of the total light is similar to most stars in this direction. This 
is also true for the magnitude of the lensed star, although its colour is unknown. Therefore 
OGLE-1999-BUL-32 is entirely consistent with being approximately at the Galactic centre. 
In the same diagram, the red clump stars around I = 15.3 and V — I = 1.8 are clearly 
visible. 

The (online) DoPhot photometry is quite noisy, because the lensed star is heavily blended 
with nearby star(s) (see §3), the fluctuations in the seeing make it difficult for DoPhot to 
disentangle the relative contributions. In fact, it is so noisy that the time-scale of this event 

t http:/ /www.astrouw.cdu.pl/" oglc/oglc2/cws/1999/bul-32.html; 
ftp://darkstar.astro.washington.edu/macho/Alert/99-BLG-22/ 
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is hard to determine with the DoPhot photometry. In contrast, the difference image analysis 
automatically subtracts out any blending. As a result, the errors are much reduced and the 
number of usable images is also increased. Both improvements are crucial for determining the 
long-duration nature of the event and, more importantly, for detecting the subtle parallax 
effect. Initially we analysed just the three season data from 1997 to 1999 available online 
(Wozniak et al. 2001). However, the parallax model predicts deviations from the standard 
model in the 2000 season and to test this we subsequently analysed the data from this season. 
Reassuringly, this confirmed the prediction of our parallax model. The four-season data from 
the difference image analysis is shown in Fig. 2§. In total, there are 246 data points in the 
light curve. In the next section, we present both the best standard and parallax model for 
this unique event. 



3 MODEL 

We first fit OGLE-1999-BUL-32 with the standard single microlens model. In this model, 
the (point) source, the lens and the observer are all assumed to move with constant spatial 

velocities. The standard light curve, A{t) is given by (e.g. Paczyhski 1986): 

7/ 2 -I- 9 i 

Ait)= T== , u(t) ee sfl + r{t)\ (1) 
uy u + 4 

where uq is the impact parameter (in units of the Einstein radius) and 

r(t) = ^, tE = f, (2) 

t E V 

with to being the time of the closest approach (maximum magnification), the Einstein ra- 
dius projected onto the observer plane, v the lens transverse velocity relative to the observer- 
source line of sight, also projected onto the observer plane, and tE the Einstein radius crossing 
time. The Einstein radius projected onto the observer plane is given by 



AGMD s x 

rE = V^r^)' (3) 

where M is the lens mass, D s the distance to the source and x = D\/D s is the ratio of the 
distance to the lens and the distance to the source. Eqs. (1-3) shows the well-known lens 
degeneracy, i.e., from a measured £e, one can not infer v, M and x uniquely even if the 
source distance is known. 

The flux difference obtained from difference image analysis can be written as 

§ The data are available at http://www.jb.man.ac.uk/" smao/bul33.dat 
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f(t) = f L [A{t) - 1] + Af, (4) 

where /l is the baseline flux of the lensed star, and Af = / — /r is the difference between the 
baseline flux (/ ) and the flux of the reference image (/«). fo includes the (unmagnified) flux 
of the lensed star and blended star(s), if present. Note that in general Af does not have to 
be zero or even positive as the reference image can be brighter than the true baseline image 
(/r > fo)- For OGLE-1999-BUL-32, the reference image flux is f R = 359.5 (Wozniak et al. 
2001). Therefore, to fit the /-band data with the standard model, we need five parameters, 
namely, /l, Af (or / ), uo,to, and £e- Best-fit parameters (and their errors) are found by 
minimizing the usual x 2 using the MINUIT program in the CERN library T 

Our attempts to fit the light curve with the standard model reveal an ambiguity. This is 
due to the degeneracy between /l, Uq and £e for a heavily blended light curve (Wozniak & 
Paczyhski 1997). In such cases, only the combinations -u ^e and fh/uo are well determined, 
but not u ,t E and f L individually. If the parameter u is left unconstrained for this event, 
then a u value close to zero is formally preferred, with x 2 = 524.4 for 241 degrees of freedom. 
However, such a perfect alignment is statistically unlikely. For illustrative purposes, in Fig. 
2, we show the best fit with uq fixed to be 0.01, which has a slightly worse x 2 — 530.7 than 
the best fit with uq left unconstrained. The fit parameters are presented in Table 1. The top 
panel in Fig. 2 shows the difference between the data points and the standard model. Clearly 
the observed light curve shows systematic deviations from the model. Quantitatively, the x 2 
per degree of freedom is about ~ 2.2, which is unacceptably large. Since the microlensing 
variability can be clearly seen over at least four years, during which time the Earth has 
moved through four orbits, it is natural to ask whether the incorporation of the parallax 
effect will remove the inconsistency. We show next that this is indeed the case. 

To account for the parallax effect, we follow the natural formalism of Gould (2000) 
and describe the lens trajectory in the ecliptic plane. This requires two further parameters, 
namely the projected Einstein radius onto the observer plane, and an angle ip in the ecliptic 
plane, which is defined as the angle between the heliocentric ecliptic x-axis and the normal 
to the trajectory (This geometry is illustrated in Fig. 5 of Soszyhski et al. 2001). Once these 
two parameters are specified, the resulting lens trajectory in the ecliptic plane completely 
determines the separation between the lens and the observer (i.e., the quantity which is 
analogous to the standard model's uq parameter from eq. 1). This allows the light curve 

^ http:/ /wwwinfo.cern.ch/asd/cernlib/ 
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Table 1. The best standard model (first row) with the impact parameter uo fixed to be 0.01 and the best parallax model 
(second row) for OGLE-1999-BUL-32. The parameters are explained in §3. 

Model t t E (day) u f L Af i/> f E (AU) X 2 /dof 

S 1365.7 ±0.08 1495.9 ± 8.7 0.01 13.48 ± 0.05 -242.5 ±0.48 550.62/241 
P 1322tg* 641±™ 0.08to,li 43.3 ±5.1 -240.9+.^ 3.38jlo;f 29. 2+Jj 257.0/239 

to be calculated; the complete prescription is given in Soszynski et al. (2001), to which we 
refer the reader for further technical details (see also Alcock et al. 1995; Dominik 1998). For 
the parameters /l, Af, Uo,to, and t E , we take the fit parameters from the standard model 
as the initial guesses, while te and ip are arbitrarily chosen for a number of combinations 
to search for any degeneracy in the parameter space. The best-fitting parameters are again 
found by minimizing the x 2 ■ Notice that in the parallax model, u and t describe the closest 
approach and the corresponding time of the lens trajectory with respect to the Sun in the 
ecliptic plane. They no longer have straightforward intuitive interpretations as analogous 
parameters in the standard model, due to geometric projections and the parallax effect. For 
example, the closest approach in the ecliptic plane is in general not the closest approach in 
the lens plane, and hence does not correspond to the peak of the light curve. 

The model parameters for the best-fitting parallax model are presented in Table 1. The 
best fit has a x 2 — 257 for 239 degrees of freedom. We found that the lens trajectory 
parameters (u and ip) are not well-specified in the ecliptic plane, very likely due to the fact 
parallax signature is only modest for OGLE-1999-BUL-32. Fortunately, the most important 
lens parameters are well constrained, in particular, we have 

f E = 29.2l|t AU, t E = 6411™ day, Af = -240.91^, f L = 43.3 ±5.1. (5) 

The Einstein radius crossing-time is about 641 days, the largest ever reported for a mi- 
cro lensing event. The projected Einstein radius on the observer plane is also very large. 
As the flux in the reference image is f R = 359.5 (Wozniak et al. 2001), one sees that the 
total baseline flux is therefore f = f R + Af = 118.6 (cf. eq. 4). The lensed star therefore 
only contributes /l//o ~ 36.5 per cent of the total baseline flux. Note the blending frac- 
tion is well constrained in the model. The baseline /-band magnitude of the lensed star is 
about 18.1 — 2.5 log(/L//o) — 19.2 mag. The lensed star was highly magnified, reaching a 
magnification of about A max ps 32 at the peak. 

The projected Einstein radius and the time-scale £e immediately allow us to derive a 
transverse velocity projected onto the observer plane 
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v = — = 79 ± 16 km s _1 . (6) 

The lens mass can be expressed as a function of the relative lens-source distance (see 
Soszyhski et al. 2001; Gould 2000), 

- ■ § (k ~ w) - ™* (s^u)" (<$L) ■ - = £ - £ ■ < 7 > 

As can be seen from this equation, the lens mass depends on the relative lens-source parallax, 
7r rc i: If the source is about 7kpc away, and the lens lies in the disc half-way between the 
observer and the source (x = 1/2), then 7r rcl ks 0.143 mas, which gives a lens mass of about 
15M ; as a comparison, for a bulge self-lensing event with D s ps 8kpc and D\ w 6kpc, then 
7r rc i 0.042 mas, which would give a lens mass of about 4.4M . However, this latter scenario 
may be less likely since the projected velocity of the lens is relatively low (see §4, Derue et 
al. 1999). In either case, the implied lens mass seems to be rather large, well beyond the 
measured mass for neutron stars (1.4M ). 



4 DISCUSSION 

We have systematically searched for parallax events in the 520 microlensing candidates 
identified using the difference image analysis (Wozniak et al. 2001). In this process, we 
have discovered an extremely long microlensing event with an Einstein radius crossing time 
t E = 641 days, the longest time-scale ever reported. The event shows a small but significant 
parallax effect caused by the motion of the Earth around the Sun. This allows one to derive 
the projected Einstein radius on the observer plane of f-& ~ 29.2AU. We emphasize that while 
some parameters are not well-constrained, the limit on te is quite robust, and it is important 
to understand why. f E is limited from below because the parallax effect is quite subtle: a 
smaller value would mean that the Earth motion makes a larger relative excursion, and 
hence the perturbation on the light curve may become too large to be compatible with 
observations, is limited from above because if it is too large, then the parallax model 
would become similar to the standard model, i.e., it will be inconsistent with the data. 
Somewhat paradoxically, had the parallax effect been smaller than observed, the projected 
Einstein radius on the observer plane would have to be even larger, implying an even larger 
lens mass. 

In this paper, we have adopted the point source approximation, ignoring the finite size 
of the lensed star. It is important to see whether the assumption is justified particularly 
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because the star was highly magnified. The finite source size effect becomes important when 
the closest approach is smaller than or comparable to the stellar radius (Gould 1994; Witt 
& Mao 1994; Nemiroff & Wickramasinghe 1994). In the source plane, the closest approach, 
d, is given by 

d = r E 1 ~ X 1 ^ 200R Q ^^, Anax>l (8) 

From the color-magnitude diagram (Fig. 1), the lensed star is likely to have a stellar radius 
no more than a few solar radii. So the closest approach is much larger than the stellar radius, 
justifying the point source approximation. 

The derived f e and £e from the fitting allow us to express the lens mass with a dependence 
on the relative lens-source parallax (see eq. 7). If we assume the source is at D s = 8kpc, 
then the lens mass only depends the parameter, x, the ratio of the distance to the lens and 
the distance to the source. The low projected velocity constrains the value of x. If the lens 
and the Sun follow the pure galactic rotation, but the source is stationary at the Galactic 
centre, then v t = 220x/(l — x) kms^ 1 . The derived transverse velocity v t ~ 79kms _1 then 
implies x ~ 0.26, which in turn gives a lens mass of 37.3M Q . In principle, a maximum 
likelihood analysis on x can be performed following Alcock et al. (1995) using the observed 
velocity information. However, such an analysis depends on the uncertain Galactic model 
(both on the mass density distribution and the kinematics of stars). We do not perform such 
a calculation here. We note, however, that our lensed star is roughly in the same direction 
as theirs and has nearly the same projected transverse velocity (75kms _1 ), so we expect to 
obtain a similar probability distribution for x; their calculation indicates a value of x which 
is slightly smaller than the naive estimate above, and this would imply a lens mass that is 
even larger. If a star with M > few M Q is still burning nuclear fuels, it will be much more 
luminous than / = 18.1. Hence, if the lens is indeed this massive, then it must be dark, and 
it follows that it is likely to be a stellar mass black hole. 

There may be a better and empirical method to test the black hole hypothesis. While the 
photometric microlensing event is now over, the astrometric microlensing signature is still 
ongoing, owing to the much slower decay of the astrometric signature as a function of the 
impact parameter (e.g. Gould 1992; Hosokawa et al. 1993; H0g, Novikov & Polnarev 1995; 
Miyamoto & Yoshi 1995; Walker 1995; Paczyhski 1998). Ignoring the Earth's motion, the 
astrometric signature follows an ellipse. The major axis and minor axis are both proportional 
to the angular Einstein radius, given by 
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fE 1 - x „„ 1 — x 8kpc 
^E = 7f ^3.7mas -f-. 9 

The predicted astrometric motion is not very well specified due to the uncertainty in the 
trajectory. Fig. 3 illustrates the prediction for the best-fitting model with x = 0.25 (6?e = 

11 mas). The origin of the astrometry is chosen to be the position of the star when the lens 
is at infinity. One sees that the astrometric motion is no longer an ellipse due to the parallax 
effect. The largest astrometric motion from the origin is 9 E /^/8 3.9 mas for this case. The 
soon-to-be installed Advance Camera for Surveys II on board HST will be an ideal instrument 
for detecting this signature. The point spread function is well sampled for this instrument, 
and it may be able to reach an astrometric accuracy as high as 0.1 mas. HST has another 
distinctive advantage over the ground based interferometers as it can resolve the blends much 
more easily. Multi-colour data from HST will also be useful for studying the colour of the 
lensed star as currently only the /-band photometry is available. However, the astrometric 
motion is quite gradual and may be confused with the star's proper motion, hence a multi- 
year monitoring project would have to be undertaken. Spectroscopic observations of the 
lensed source will be useful to put further constraints on the lensing kinematics involved. 
The lens may also be accreting interstellar gas, and could be luminous in the X-ray. It would 
be very interesting to obtain a deep image using sensitive X-ray satellites such as XMM- 
Newton and Chandra. It will also be interesting to see whether the source is luminous in 
the radio. Radio observations have distinct advantages, as it is not affected by dust, and 
VLBI observations could reach ~ mas astrometry. The inset in Fig. 3 shows the position of 
the lens relative to the source centroid, which already reaches tens of mas for our example. 
Such a shift, if detected, will be a dramatic confirmation of our model. 
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V-I 



Figure 1. The colour-magnitude diagram for the 3.5' X 8' stellar field around OGLE-1999-BUL-32. The short-dashed line 
indicates the /-band baseline magnitude for the lensed star, while the two long-dashed lines indicate the magnitude and the 
colour of the total light from the lensed star plus nearby blend(s). 
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Figure 2. The /-band light curve for OGLE-1999-BUL-32 from difference image analysis. The solid and dotted lines are for the 
standard and parallax fits, respectively. The short-dashed line shows the baseline flux of the lcnsed star while the long-dashed 
line shows the total baseline flux of the lensed star and nearby blend(s). The approximate /-band magnitudes are indicated for 
these two baselines, together with the peak /-band magnitude. The inset shows the the light curve close to the peak. The top 
panel shows the residual flux (the observed data points subtracted by the standard model). Clearly the standard model shows 
systematic discrepancies. The curved solid line shows the prediction of the parallax model. 
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Figure 3. The predicted astrometric motion of the light centroid in the ecliptic plane. The size of the motion is proportional 
to the angular Einstein radius, which we have taken to be 11 mas (see eq. 9). The solid dots indicate the centroid positions 
every six months on and after July 31, 2001 while the open dots indicate the positions every six months before the date. The 
plus sign indicates the source position when unlensed. Notice that the scales on the two axes are different. The inset shows the 
lens position relative to the centroid of the source in mas. The open and filled circles have the same meaning as in the main 
panel. 
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